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We show that inflationary cosmology may be used to test the statistical 
predictions of quantum theory at very short distances and at very early times. 
Hidden- variables theories, such as the pilot-wave theory of de Broglie and Bohm, 
allow the existence of vacuum states with non-standard field fluctuations ('quan- 
tum nonequilibrium'). We show that inflationary expansion can transfer mi- 
croscopic nonequilibrium to macroscopic scales, resulting in anomalous power 
spectra for the cosmic microwave background. The conclusions depend only 
weakly on the details of the de Broglie-Bohm dynamics. We discuss, in partic- 
ular, the nonequilibrium breaking of scale invariance for the primordial (scalar) 
power spectrum. We also show how nonequilibrium can generate primordial 
perturbations with non-random phases and inter-mode correlations (primordial 
non-Gaussianity). We address the possibility of a low-power anomaly at large 
angular scales, and show how it might arise from a nonequilibrium suppression 
of quantum noise. Recent observations are used to set an approximate bound 
on violations of quantum theory in the early universe. 
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1 Introduction 



According to inflationary cosmology [T], the early universe underwent a period 
of exponential expansion, during which microscopic quantum fluctuations were 
stretched to macroscopic scales. The resulting (classical) primordial pertur- 
bations seem to be of the form required to explain the observed temperature 
anisotropy in the cosmic microwave background (CMB), and are widely be- 
lieved to have seeded the formation of large-scale structure generally. In this 
scenario, precision CMB measurements today can provide information about 
— and tests of — microscopic physics in the very early universe. For this rea- 
son, many workers have turned to inflationary CMB predictions in the hope 
that these will provide a 'cosmic microscope' with which to probe high-energy 
physics at very short distances and at very early times. However, if the pri- 
mordial perturbations do indeed have a quantum origin, then inflationary CMB 
predictions will also be sensitive to the structure of quantum theory itself, as well 
as to that of high-energy physics. Therefore, inflationary cosmology and CMB 
measurements may equally be used to probe possible deformations of quantum 
theory at very short distances and at very early times. 

In a typical inflationary scenario, at very early times the cosmological scale 
factor a{t) undergoes a period of approximately exponential growth, a oc e^* 
with H « const.. During inflation, field perturbation modes have physical 
wavelengths Aphys — a{t)X (x e^*. (As usual, A — 27r/fc is the wavelength 
today — the 'comoving wavelength' — and we set the scale factor today to be 
uq — 1.) A mode 'exits' the Hubble radius H^^ when Aphys ^ H~^, at a time 
texit = icxit(fc) (which can be defined by 27ra(toxit)/^ ~ or by a(icxit)/^ ~ 
H~^). Soon after ioxit(^), the perturbation 'freezes' and becomes part of the 
primordial spectrum. After inflation ends, physical wavelengths Aphys oc a grow 
more slowly than the Hubble radius H^^ = a/a cx t (where a oc t-^^^ or for 
radiation-dominated or matter-dominated expansion respectively). Mode 're- 
entry' occurs at a time icntcr(fc) when Aphys ^ H~^, after which the (formerly- 
frozen) perturbations begin to grow, eventually giving rise to anisotropics in the 
CMB and to large-scale structure l2j . 

While there are many uncertainties surrounding the details of inflationary 
cosmology, there is a broad consensus that the formation of (frozen) primordial 
perturbations takes place when the corresponding physical wavelengths Aphys ^ 
H^^ are truly microscopic. Further, because of the huge expansion during 
the inflationary phase, the relevant modes will have had very short physical 
wavelengths, Aphys << at the onset of inflation (where the shorter the 

wavelength, the later the time at which the mode exits the Hubble radius during 
the inflationary phase.) Indeed, it appears that even modes with initial Aphys ^ 
^P, where Ip « lO^'^'^ cm is the Planck length, may contribute to the primordial 
spectrum [3]. Clearly, if inflation did indeed occur, then precision measurements 
of the CMB (and of large-scale structure generally) can probe physics at very 
early times and at very short distances (possibly even at distances < Ip, to the 
extent that this might be meaningful). 

A number of possible deformations of high-energy physics have been consid- 
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ered in an inflationary context. These include: (a) modified dispersion relations 
(which may be introduced ad hoc [3l|4j, or which may be motivated by quantum- 
gravitational deformations of Lorentz invariance [5] or by quantum cosmology 
[B]); (b) an ultra-violet cutoff coming from a fundamental length associated 
with deformed uncertainty relations (possibly associated with quantum gravity 
or string theory) [7j; (c) short-distance non-commutative geometry [5]. Some 
authors consider that changes in physics at very high energies may have an ef- 
fective description in terms of different choices of quantum vacuum [3] (for a 
review, see ref. [TO])- Excited, non- vacuum states have also been considered [TT]. 
However, while it is generally agreed that inflationary primordial perturbations 
have a quantum origin, effects on the CMB arising from possible deformations 
of quantum theory itself are not usually considered. (By 'quantum theory' we 
mean, essentially, the representation of physical states in Hilbert space, with 
unitary evolution, and with probabilities given by the Born rule.) An exception 
is Perez et al. [12], who discuss how predictions for the CMB could be affected 
by a hypothetical dynamical collapse of the wave function, a proposal that is 
motivated by the quantum measurement problem (which seems especially severe 
in a cosmological setting). 

Despite the widespread reluctance to consider deformations of quantum the- 
ory itself, there is in fact no good scientific reason for believing that the structure 
of standard quantum theory is 'final', or that the predictions of quantum theory 
will continue to hold under all conditions. The following arguments are often 
presented as evidence for the finality of quantum theory: that it provides a 
universal framework applicable to all systems independently of their composi- 
tion (electrons, fields, atoms etc.); that it is based on simple, elegant axioms; 
that it provides the basis for powerful new technologies; and, of course, that in 
all cases so far it agrees with experiment. However, arguments similar to these 
could have been made in the eighteenth and nineteenth centuries concerning the 
status of Newtonian mechanics: at that time, Newtonian mechanics seemed to 
provide a universal framework applicable to all systems independently of their 
composition (rocks, fluids, planets etc.); it was based on simple, elegant axioms 
(Newton's three laws of motion); it provided the basis for powerful new tech- 
nologies; and it agreed with all experiments performed to date. And yet, we now 
know that Newtonian mechanics is in fact merely approximate and emergent, 
arising from a classical and low-energy limit of relativistic quantum field theory. 
Of course, that Newtonian mechanics proved to be approximate and emergent 
does not imply that quantum theory will necessarily turn out likewise. However, 
the case of Newtonian mechanics does suggest that the above (frequently-cited) 
arguments for the finality of a physical theory are not reliable. 

The ultimate test of the domain of validity of a scientific theory is, of course, 
experiment. No matter how well a theory has been tested in the past, it will 
always be subject to possible modification in the future, in hitherto untested 
regimes. Therefore, in order to expand our knowledge of the domain of validity 
of any given theory, it is necessary to subject it to ever more stringent tests in 
ever more extreme conditions. To accomplish this, it is helpful to have a 'foil' 
against which to test the theory in question — that is, to have a model reducing 
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to the given theory only in some hmit. 

In the case of quantum theory, a number of alternatives or foils might be con- 
sidered. Models with a nonlinear evolution or with a dynamical collapse of the 
wave function have, for example, been subjected to considerable experimental 
scrutiny. In this paper, we focus on a different possibility: that of nonequilib- 
rium hidden variables [Il[Il[I3[Tni[I71[TSl[Tl[inilIIl[^ll3[Mlll5]. 

A deterministic hidden-variables theory, such as the pilot-wave theory of 
de Broglie [17| and Bohm [5H], agrees with quantum theory only in the 
limit in which the hidden parameters have a particular 'quantum equilibrium' 
distribution [T31 [TSl [THl [H] • A foil against which to test quantum theory may 
then be obtained from such a theory by allowing the hidden variables to have a 
non-standard or 'nonequilibrium' distribution, resulting in statistical predictions 
that deviate from those of quantum theory [25]. Such possible corrections to 
quantum theory will be explored here, in the context of inflationary cosmology, 
where it will be shown how CMB observations may be used to set bounds on 
the presence of quantum nonequilibrium at very short distances and very early 
times. 

In section 2, we review the notion of quantum nonequilibrium, in de Broglie- 
Bohm theory and in general (deterministic) hidden-variables theories, and we 
provide motivation for why quantum nonequilibrium might exist in the very 
early universe. In section 3, we develop pilot-wave field theory on an expanding 
space, and we write down equations for the time evolution of arbitrary (nonequi- 
librium) distributions in an expanding universe. In section 4, we discuss two 
scenarios whereby quantum nonequilibrium could exist during inflation: first, 
nonequilibrium for large- wavelength modes could survive from a pre-inflationary 
era, since under the right conditions relaxation can be suppressed at large wave- 
lengths on an expanding space; second, nonequilibrium might be generated by 
novel gravitational processes at the Planck scale. In section 5, we review the 
standard theory of CMB temperature anisotropies, their explanation in terms 
of primordial curvature perturbations, and the production of the latter by infla- 
ton fluctuations during inflation. In section 6, we calculate the time evolution 
of quantum nonequilibrium in the Bunch-Davies vacuum on de Sitter space, 
and we show that the width Dk {t) of the nonequilibrium distribution for each 
mode of wave number k remains in a fixed ratio \/^{k) = Dk{t)/ Ak{t) with 
the equilibrium (quantum) width Ak(t). In section 7, we show how the power 
spectrum for the primordial curvature perturbations is corrected by the factor 
^(fc). Some general remarks are made in section 8, concerning the transfer of 
microscopic nonequilibrium to cosmological scales, the effective quantum mea- 
surement of the inflaton field during the 'quantum-to-classical' transition, and 
the weak dependence of our results on the details of pilot-wave dynamics. In 
section 9, we use current CMB data to derive an approximate bound on quan- 
tum nonequilibrium during inflation; specifically, under certain assumptions, we 
show that the hidden- variable relative entropy Shv{k) (which measures the dif- 
ference between nonequilibrium and quantum probabilities for a mode of wave 
number k) satisfies the approximate bound |S'hv(fc)| 10^^ for values of k 
close to fco = 0.002 Mpc~^. In section 10, we consider the possibility of a low- 
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power anomaly at large angular scales, and we discuss how it might arise from 
a nonequilibrium suppression of quantum noise (^(fc) < 1) in certain regions 
of A:-space. In section 11, we show how nonequilibrium can generate primor- 
dial perturbations with non-random phases and inter-mode correlations. Our 
conclusions are given in section 12. 

2 Quantum Equilibrium and Quantum Nonequi- 
librium 

The notion of quantum nonequilibrium was first discussed in detail in terms of 
de Broglie-Bohm theory [13l [141 E] i and was later generalised to include all 
(deterministic) hidden- variables theories [18l [191 ISH [24j ■ 

Consider, for example, the very simple case of de Broglie-Bohm theory ap- 
plied to a single nonrelativistic particle with mass m and no spin. The wave 
function -ip = 1-01 e^^ (with units ft = 1) acts as a 'pilot wave' that deter- 
mines the velocity of the particle according to de Broglie's guidance equation 
dx/dt — (l/m)VS' (or dx/dt = (1/m) Im (Vi/'/^/')) — an equation that deter- 
mines the trajectory x(f) of the particle, given the initial position x(0) (assuming 
that ■0 = V'(^j t) is known for all x and t, by solving the Schrodinger equation 
with a given initial wave function V'(x,0)). Let "0 propagate in free space, then 
strike a screen with two slits, and finally strike a backstop where the particle 
is detected. The pilot wave undergoes interference upon traversing the screen. 
The location x(i) of the particle at any time t is determined (in principle) by 
the initial value x(0); in particular, where the particle lands on the backstop 
is determined by x(0). Because the velocity field (l/m)VS' is equal to the 
usual quantum probability current j divided by the usual quantum probabil- 
ity density iV'l^i it follows trivially that an initial ensemble of particles guided 
by the same pilot wave ip and with positions x(0) distributed according to the 
equilibrium rule p{x, 0) = |-0(x, 0)| will evolve into an equilibrium distribution 
pipi.t) — ['(/'(x, t)| at later times, resulting in the usual quantum distribution 
of particles at the backstop (showing the usual interference pattern). On the 
other hand, it is easy to see that in general an initial 'nonequilibrium' ensem- 
ble with distribution p(x, 0) 7^ |0(x, 0)|^ results in a non-quantum distribution 
p(x, i) 7^ |-0(x, t)| at the backstop. (For example, in the absence of a rapid di- 
vergence of neighbouring trajectories, if p(x, 0) is concentrated around a single 
initial point x(0) then p(x, t) will be concentrated around a single trajectory 
x(t), and the usual interference pattern will be replaced by a single localised 
spot.) 

The pilot-wave theory of a many-body system was first proposed by de 
Broglie at the 1927 Solvay conference [26l [27]. For a system of n (nonrela- 
tivistic) particles with positions Xi(t) and masses rrii, de Broglie's law of motion 



takes the form 



dxi 1 
—— ~ — Im 

dt rrii 



Tp mi 



(1) 
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where V' = x„, i) is the many-body wave function. De Broghe regarded 

([1]) as expressing a unification of the principles of Maupertuis and Fermat, re- 
sulting in a new form of dynamics based on velocities ^7} . 

Writing the total configuration as q = (xi, X2, x„), it is again readily 
shown that for an ensemble of systems guided by the same wave ip and with 
configurations distributed according to p{q,0) = |V'('Z, 0)|^, the distribution of 
configurations at later times will be p{q,t) = \ip{q,t)\ . 

As shown in detail by Bohm in 1952 [5S], the above 'de Broglian' dynam- 
ics may be applied to the process of quantum measurement itself, by treating 
the system being measured together with the measuring apparatus as a single 
many-body system of n particles. The total configuration q = (xi,X2, ....,x„) 
then defines the 'pointer position' of the apparatus, as well as defining the con- 
figuration of the measured system. For each run of a quantum experiment, the 
evolution is deterministic: the initial conditions <7(0), ■(/'(q, 0) determine the final 
conditions q(t), '4>{q,t). Over an ensemble of initial configurations q{0) guided 
by the same wave function -0, if we assume the initial quantum equilibrium con- 
dition p{q,0) = \il'{q,0)\^, then the statistical distribution of pointer positions 
at later times will agree with the predictions of quantum theory. 

Schematically, during a standard quantum measurement, the initial packet 
ipiq, 0) on configuration space evolves into a superposition tp{q, t) = J2n 
of terms ipnilT^) that separate with respect to the pointer degrees of freedom 
(that is, distinct ip^{q,t) have negligible overlap with respect to the pointer de- 
grees of freedom). The final configuration q{t) can then be in (the support of) 
only one 'branch' of the superposition, say ip^ {q, t) . For an initial equilibrium 
ensemble, it is readily shown that this occurs with probability \ci\ , in accor- 
dance with the Born rule. Further, inspection of de Broglie's velocity law ([1]) 
shows that the motion of q{t) will then be affected by ip^{q,t) alone, resulting 
in an effective 'reduction' of the wave function. 

As in the simple example of a single particle undergoing interference, for a 
general quantum measurement the distribution of outcomes depends crucially 
on the assumed initial distribution p{q,0) of initial configurations q{0). For 
a nonequilibrium ensemble, p{q,0) 7^ |';/'(9, 0)|^, the distribution of quantum 
measurement outcomes will generally disagree with the predictions of quantum 
theory (assuming that relaxation to equilibrium has not taken place in the 
meantime — see below). 

De Broglie's dynamics may be readily applied to fields, where (say for a 
scalar field cj}) the motion of the field configuration q{t) ~ t) is determined 
by the Schrodinger wave functional 5'[0(x),t]. Indeed, for any system with 
configuration q, as long as the Schrodinger equation idip/dt = Hip for ipi{q,t) 
has an associated locally-conserved current j — j [ip] — j{q,t) in configuration 
space, obeying a continuity equation 



V • J - 



dt 

(with V = d/dq), one may define a de Broglian or pilot- wave dynamics for the 
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system, by introducing the configuration-space velocity field 

dq ^ _J_ 
dt ~ ■ 

(Note that in this dynamics, ip is viewed as a physical field in configuration space, 
guiding the motion of an individual system, with no a priori connection with 
probabilities.) For an ensemble of systems, each with the same wave function 
ip{q,t), we may consider an arbitrary initial distribution p{q,0) ^ \ip{q,0)\ , 
whose time evolution p{q, t) is determined by the de Broglian velocity field q in 
accordance with the continuity equation 

| + V.(p^)^0. 

Because obeys the same equation, an initial distribution 0) — 0)|^ 
evolves into p{q,t) = \il:{q,t)\^ . This is the state of quantum equilibrium, but 
the theory clearly allows one (in principle) to consider any initial distribution 
— just as classical mechanics allows one to consider any initial distribution 
departing from thermal equilibrium. 

It is worth emphasising that pilot-wave theory is a radically new form of 
dynamics, very different from classical (Newtonian or Hamiltonian) mechanics. 
This was in fact de Broglie's original point of view, but it was unfortunately 
obscured by Bohm's pseudo-Newtonian reformulation of the theory in terms of 
a law for acceleration (involving a 'quantum potential') [27j . 

Pilot-wave dynamics is grounded in configuration space, where ip propagates. 
While the dynamics is local in configuration space, it is highly nonlocal when 
projected down to 3-space (as required by Bell's theorem). For example, if a 
particle with position xi is entangled with a particle with position X2, then 
the velocity xi depends instantaneously on X2 (no matter how remote X2 may 
be from xi), and changing the local Hamiltonian at X2 is found to have an 
instantaneous effect on the distant velocity xi . Such nonlocal effects are erased 
upon averaging over an equilibrium ensemble p = but in nonequilibrium, 
p ^ \ip\ , there are (in general) nonlocal signals at the statistical level p!4l [T5] . 
suggesting the existence of an underlying preferred foliation of spacetime [29] . 

Pilot-wave dynamics — as originally formulated by de Broglie — is also 
first order in time in configuration space (rather than in phase space): the 
fundamental law of motion determines velocities, not accelerations. This last 
feature has important implications for the associated kinematics: for particles, 
the natural state of motion is rest (instead of uniform motion in a straight line) , 
and there is indeed a natural preferred foliation of spacetime with a fundamental 
time parameter t (consistent with the fundamental nonlocality of the theory) 

ISQi. 

Quantum nonequilibrium may be considered, not only in pilot-wave the- 
ory, but also in any deterministic hidden- variables theory [181 IHl IMj ■ For in 
any such theory, given macroscopic experimental settings M, there is a mapping 
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Lo ~ uj{M, A) from initial hidden variables A to final outcomes w of quantum mea- 
surements. There is also a 'quantum equilibrium' probability measure Pqt('^)> 
defined on the set of hidden variables, that yields quantum probabilities Pqt(w) 
for the outcomes. (In the case of pilot-wave theory, Pqx(A) is given by p = |'0|^-) 
Once such a theory has been constructed, one may consider arbitrary 'nonequi- 
librium' probability measures p(A) ^ PqtI''^), resulting in outcome probabilities 
P{i-u) 7^ Pqt{uj) that depart from the predictions of quantum theory. 

In this paper we shall be studying quantum nonequilibrium in the context of 
inflationary cosmology, using the pilot-wave theory of fields as a concrete exam- 
ple. However, we emphasise that similar studies could be made in any determin- 
istic hidden- variables theory, simply by making the replacement /9qx(A) — *■ p{X)- 

At present, pilot- wave theory is the only deterministic hidden- variables the- 
ory of broad scope that we possess, though some attempts have been made to 
construct alternative theories. For example, in the 1980s, Smolin attempted to 
construct a deterministic hidden-variables theory of an A'^-body system, based 
on the classical Hamiltonian dynamics of a certain N x N matrix Mij (t) , whose 
eigenvalues correspond to particle positions and whose off-diagonal elements 
correspond to nonlocal hidden variables associated with pairs of particles [31| . 
Adopting a classical action principle for the (deterministic) dynamics of the 
matrix, Smolin made a number of assumptions, including a statistical assump- 
tion to the effect that the coarse-grained evolution of the off-diagonal terms 
amounts to a Brownian motion. In the limit of a large number A^ of particles 
with masses m^, it was shown from these assumptions that the particle positions 
also undergo a Brownian motion, that the ith particle current velocity (the 
average of the mean forward and backward velocities) is given by a gradient, 
Vi = ^iS/rrii, where 5* is a function on configuration space, and that the com- 
plex function ip = ^/pe*'^ (where p is the particle probability distribution on 
configuration space) satisfies the Schrodinger equation for a many-body non- 
relativistic system. Smolin's strategy was to show that his assumptions led, in 
the limit of large A^, to the basic postulates of Nelson's stochastic mechanics 
[32]. As was already known, in Nelson's theory — which is based on a form 
of Brownian motion subject to special conditions, including the condition that 
Vi = ViS/rrii for some function S — the derived quantity ip = y/pe^^' indeed 
satisfies the Schrodinger equation. 

More recently, a model similar to the above (though based on the bosonic 
part of the classical matrix models used in string and M theory) was again inves- 
tigated by Smolin, with similar assumptions and results [33j . In ref. [31j , it had 
also been suggested that one might consider a model in which the off-diagonal 
matrix elements of Mij{t) are constant, with fluctuations in a local system aris- 
ing from the nonlocal transmission of fluctuations from other particles in remote 
regions of space. This last model has recently been recast in terms of the dy- 
namics of a graph with A^ nodes [51]: assuming that the edges of the graph do 
not evolve in time, the corresponding adjacency matrix is constant, and is taken 
to be the off-diagonal part of matrices Mij(t). Again, as in Smolin's original 
model, assumptions are made so as to arrive at Nelson's stochastic mechanics 
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in some approximation. 

However, while it is often claimed that Nelson's theory is empirically equiv- 
alent to quantum theory, unfortunately, as shown by Wallstrom [3S], the two 
theories are in fact not equivalent, because Nelson's function S does not have 
the specific multivalued structure required for the phase of a single- valued (and 
continuous) complex field ip. The Schrodinger equation is indeed derived, but 
only for the exceptional set of wave functions with no nodes, for which the circu- 
lation of Vi5' around all closed curves vanishes. Since almost all wave functions 
have nodal points (where ip = 0), quantum theory cannot be derived from Nel- 
son's theory, or from any model that leads to Nelson's theory. (Note that there 
is no such problem in pilot- wave theory, where ip is regarded as a basic entity.) 

Thus, as they stand, the deterministic models of refs. [STJ [33l [34] seem to 
yield derivations of Nelsonian mechanics, but not of quantum mechanics. Some 
basic element is missing. One must somehow ensure that the circulation of V^S 
around nodes of p can be non-zero but always restricted to integer multiples of 
27r. (And if one wishes to derive the wave function, then of course one cannot 
simply assume at the outset that S is the phase of a complex-valued field.) 
Still, if some way were found to solve Wallstrom's phase problem, then such 
derivations of Nelsonian mechanics as an average over a certain statistical state 
could again be generalised to arbitrary statistical states, yielding nonequilibrium 
departures from quantum theory in the sense considered heren 

As another example, Adler [35] has constructed what appears to be a deter- 
ministic hidden-variables theory, in which the parameters A are matrices with 
Grassmann (even and odd) valued elements, obeying a generalised form of clas- 
sical Hamiltonian dynamics. The state of thermal equilibrium, defined in the 
usual way on phase space, is argued to lead (after some approximations) to 
a quantum-like phenomenology with a dynamical wave function collapse. The 
precise nature of Adler's theory seems to require further elucidation; but if it is 
indeed a hidden- variables theory in the sense meant here, then thermal nonequi- 
librium in Adler's theory should again correspond to quantum nonequilibrium. 

In the author's view, because Hamiltonian dynamics is of second order in 
configuration space, it is not a natural framework for nonlocal theories with 
a preferred state of rest or preferred slicing of spacetime — unlike pilot-wave 
dynamics, which is first order in configuration space, and which therefore (as 
we have mentioned) provides a natural setting for such theories |30) . But even 
so, the above alternative theories based on Hamiltonian dynamics do illustrate 
that the idea of quantum nonequilibrium is a general one. 

For the purposes of this paper, it suffices that there exists at least one 
model of quantum nonequilibrium, based on pilot-wave dynamics, that may 

■^Smolin ,36^ has attempted to solve Wallstrom's phase problem by allowing discontinuous 
wave functions. However (even leaving aside the resulting divergences for expectation values 
of quantum observables such as kinetic energy), Smolin applies his prescription only to the 
case of a particle moving on a circle, which is too simple to capture the nature of the problem 
raised by Wallstrom. In higher dimensions — for example even in two dimensions, and with 
just one node — allowing discontinuous wave functions results in an ill-defined (one-to-many) 
mapping from Nelsonian states to quantum states. For a full discussion, see ref. |37| . 



9 



serve as a foil against which to test quantum theory. To be able to provide 
quantitative bounds on violations of quantum theory in the early universe is 
motivation enough to consider models with quantum nonequilibrium. Even 
so, before proceeding, let us briefly provide some further motivation for why 
quantum nonequilibrium might exist at very early times. 

First, it has been shown that in pilot-wave theory the equilibrium state p = 
I'i/'l^ may be understood as arising from a process of relaxation that is analogous 
to classical thermal relaxation, where the former is defined on configuration 
space rather than on phase space. The difference between p and \ip\^ may be 
quantified by the i/-function 

j dgpln(p/|^|2) (2) 

(equal to minus the relative entropy of p with respect to li/"]^), which obeys a 
coarse-graining i/-theorem analogous to the classical one, and where the mini- 
mum H = corresponds to p = \tp\^ |131ll5l[r7] . Further, numerical simulations 
for simple two-dimensional systems [53J[3n] show a remarkably efficient approach 
to equilibrium, with an approximately exponential decay of the coarse-grained 
77-function, H{t) — > 0, and a corresponding coarse-grained relaxation p — > 1-01^ 
(assuming appropriate initial conditions for p and ip)^ Because all the systems 
we have access to (such as hydrogen atoms in the laboratory) have a long and 
violent astrophysical history, we would then expect to see quantum equilibrium 
in these systems. While it is logically possible, of course, that the universe was 
simply born in a state of quantum equilibrium, it seems more natural to con- 
sider that the equilibrium we see today arose from relaxation processes in the 
remote past [16| 117) , in which case the very early universe is the natural place 
to look for nonequilibrium phenomena. 

Second, an appealing feature of this picture concerns the status of locality in 
physics. It may be shown that quantum nonequilibrium for entangled systems 
leads to nonlocal signals at the statistical level, in pilot- wave theory (as already 
mentioned) and indeed in any deterministic hidden-variables theory; while in 
equilibrium, the underlying nonlocal effects cancel out at the statistical level 
[141 US [HI US] ■ Locality is therefore a contingency (or emergent feature) of 
the equilibrium state. Similarly, standard uncertainty-principle limitations on 
measurements are also contingencies of equilibrium jT4l [151 UHl HI] • These re- 
sults provide an explanation for the otherwise mysterious 'conspiracy' in the 
foundations of current physics, according to which (roughly speaking) quantum 
noise and the uncertainty principle prevent us from using quantum nonlocality 
for practical nonlocal signalling. From the above perspective, this 'conspiracy' 
is not part of the laws of physics, but merely a contingent feature of the equilib- 
rium state (much as the inability to convert heat into work, in a state of global 
thermal equilibrium, is not a law of physics but a contingency of the state). On 

^The understanding of relaxation in pilot-wave theory is subject to the usual caveats — 
familiar from classical statistical mechanics — associated with initial conditions and time 
reversal. For detailed discussions of this point, see refs. 1151 1161 [TTl I23| . 
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this view, quantum physics is merely the effective description of a particular 
state — just as, for example, the standard model of particle physics is merely 
the effective description of (perturbations around) a particular vacuum state 
(arising from spontaneous symmetry breaking). If one takes this view seriously, 
it suggests that nonequilibrium phenomena should exist somewhere (or some 
time) in our universe. And again, the early universe seems the natural place to 
look. 

Quantum nonequilibrium at very early times may also be motivated by the 
cosmological horizon problem, which may be avoided by the explicit nonlocality 
associated with nonequihbrium [TH [TSl HH US] — see section 4.1. 

Finally, if one takes de Broglie-Bohm theory seriously, one should take the 
possibility of nonequilibrium seriously as well, since it is only in nonequilibrium 
that the underlying details of the theory become visible (via measurements 
more accurate than those allowed by quantum theory [201 El]). If instead the 
universe is always and everywhere in quantum equilibrium, the details of de 
Broglie-Bohm trajectories will be forever shielded from experimental test, and 
de Broglie-Bohm theory itself would be unacceptable as a scientific theory. 

For the above reasons, then, we are led to consider the hypothesis of quan- 
tum nonequilibrium at or close to the big bang [T31 [HI [13 [HI [13 [IS] . It is 
the purpose of this paper to show that inflationary cosmology provides a means 
of testing this hypothesis, through precision measurements of the cosmic mi- 
crowave background. 



3 Pilot- Wave Field Theory on Expanding Space 



where again a{t) is the scale factor, with Hubble parameter H = a/a. As is 
customary, we take ao = 1 today (at time to), so that |dx| is a comoving distance 
(or proper distance today). 

A free (minimally-coupled) massless scalar field <j> has a Lagrangian density 



For simplicity we restrict ourselves to a fiat metric. 



(3) 




with an action J dt J d^x C (where x are comoving coordinates). This implies 
a canonical momentum density tt = dC/dcf) = a^cj) and a Hamiltonian density 



The equations of motion = SH/Sir, n 
the classical wave equation 



dH/S(l) (with H = Jd^x H) lead to 



+ 



a a^ 



(6) 
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Pilot-wave field theory is defined in terms of the functional Schrodinger pic- 
ture, with a preferred fohation of spacetime [l5 l fT6 l [28 l l40 l[4Tl[42l[43ll44l l45]. 

For an expanding universe with metric ([3]), containing a scalar field (f) with 
Hamiltonian density a general wave functional ^'[0, t] = (0(x)|4'(i)) (where 
|(/)(x)) is a field eigenstate) satisfies the functional Schrodinger equatioro 



dt J \ 2a? Scf)'^ 2 

(with the usual realisations — s- 0, tt — > —iS/Scf)). This implies the continuity 
equation 

(where 5* = |4'| e*"^), from which one may identify the de Broglie velocity 

dt a3 6<P ^ ' 

for an individual field configuration. Here, again, ^ is interpreted as a physical 
field in configuration space, guiding the evolution of an individual field 0(x, t) in 
3-space. (Note that S is defined only locally, as 5 = Imln^'. One may equally 
write ([U as ^ = ^Im-^|^, without mentioning S.) 

A similar construction may be given in any globally-hyperbolic spacetime, by 
choosing a preferred foliation [27 . Thus there is no need for spatial homogeneity. 

Over an ensemble of field configurations guided by the same pilot wave VP, 
there will be some (in principle arbitrary) initial distribution P[(j),ti], whose 
time evolution t] will be determined by 

^d^.l(pV4)-0. (10) 



dt J (50 V 

UP[(j),t.i] = then P[(j>,t] — |\E'[(/), for alH, and empirical agreement 

is obtained with standard quantum field theory [28l HH US [43l [Ml |45] . On the 
other hand, for an initial nonequilibrium distribution P[(j>,ti] ^ |\E'[0, for 
as long as P remains in nonequilibrium the predicted statistics will generally 
differ from those of quantum field theory. In any case, whatever form P may 
take (equilibrium or nonequilibrium), its time evolution will be given by pop . 

It will prove convenient to rewrite the dynamics in Fourier space. Expressing 
(^(x) in terms of its Fourier components 



1 



"^•^ (2^)3/2 

and writing 



d^x <hMe-'^-'- 



(27r)3/ 

*As usual in this context, some sort of regularisation is implicitly assumed. 
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for real q^r {r = 1, 2), where V is a box normalisation volume, the Lagrangian 
L = J (Px£ becomes 



2 

kr 

(For V ^ oo, ^Y.k^ (S^F -f ^^'^ ^^kk-^ {2TT)^S^(k - k). The reality of 
(j) requires (j)^ = <p_]^ or g^i = q-ki, qk2 = —Q-k2, so that a sum over physical 
degrees of freedom should be restricted to half the values of k, for example 
kz > 0.) Introducing the canonical momenta 

dL 



Tl'kr = TT. = a 9kr , 

oqkr 



the Hamiltonian becomes 



1 

,2a3 ' 2' 

kr 

The Schrodinger equation for ^ = '^[qkr, t] is then 

^2 



^^ = E(-^ai; + 2"'^-j*' ^''^ 

which implies the continuity equation 



(13) 



dt ^ dqur \ dqur J 

kr 

and the de Broglie velocities 

dqkr 1 dS 



dt dqkr 



(again with ^ = \^\e^^). The time evolution of an arbitrary distribution 
P[qkr,t\ will then be given by 

For product states 



*[9kr,i] =n^kr(9kr,t) (15) 
kr 

(such as the Bunch-Davies vacuum during inflation), the wave function tp^^^ for 
a single mode kr satisfies 



.ay;,,, ^ / 1 9^ 1 

dt V 2a^dql^^2 



i-^ ={-^^ + T^ak^lir ^kr • (16) 
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Writing V'kr = IV'krI s'"'"' (where S = J2kr *kr), the de Broghe velocity for 
is then 

dqkr 1 9skr 



dt dqi^r 

If the initial distribution P[qkr,ti\ also takes the product form 



(17) 



P[qkr, k] = Yl Pkrl'Zkr, U) , (18) 



kr 



then the time evolution of />kr(9kr7^) wih be given by 



''^^^ ■ ' ^P.AP^]-0. (19) 



Note that the factorisability condition (|T8|) for the probability distribution 
P is logically independent of the factorisability condition (fT5|) for the pilot wave 
5". Thus, even for a vacuum state, in nonequilibrium it is still possible to have 
inter-mode correlations. For simplicity, in section 6, we shall restrict ourselves to 
the case of uncorrelated nonequilibrium modes. The correlated case is discussed 
in section 11. 



4 Quantum Nonequilibrium in the Very Early 
Universe 

In this paper, the focus is on setting experimental bounds on possible violations 
of quantum theory during inflation. Before proceeding with this, however, let 
us indicate how one might (in future work) be able to predict details of such 
violations. The scenarios sketched in this section also serve to give a preliminary 
idea of the kinds of violations one might expect to find. 

4.1 Relic Nonequilibrium from a Pre-Inflationary Era 

One reason to expect early nonequilibrium to exist is that, as sketched in sec- 
tion 2, according to de Broglie-Bohm theory ordinary matter corresponds to a 
'quantum equilibrium phase', and it is natural to suppose that this equilibrium 
state emerged from the violence of the big bang. 

Another reason is that nonequilibrium at very early times would unleash the 
nonlocality inherent in all hidden- variables theories, thereby evading the horizon 
problem associated with an early Friedmann expansion (if there was one). For 
a oc i^/^ the horizon distance is (with c = 1) 

hit) = ait) f ^ = 2t , 
Jo o-it) 

and for any two comoving points separated by a coordinate distance |Ax| we 
have Ihit) « ait) |Ax| for sufficiently small t. On this basis it has been widely 
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argued that early homogeneity — over seemingly causally-disconnected domains 
— is unnatural and puzzlingjj As we have mentioned, the hypothesis of quan- 
tum nonequilibrium at the big bang was originally introduced partly to solve 
this problem [TH [1^1 HH HH] • For the above scalar field, for example, a generic 
wave functional will be entangled across space, so that the field velocity 

(90(x,i) _ 1 dS[(j),t] 
dt ~ ^ (5(/)(x) 

at a point x will depend on instantaneous values of the field at remote points 
x' 7^ X, and in nonequilibrium this nonlocal dependence will not be hidden by 
statistical noise (as it is in quantum theory). Of course the horizon problem 
was also one of the historical motivations for introducing inflation: the period 
of exponential expansion ensures that our observable region originates from 
within a single causal patch [46]. However, even in an inflationary context, it 
appears that some models require homogeneity as an initial condition in order 
for inflation to begin [47]. Therefore, it is possible that consideration of a pre- 
inflationary era will revive the horizon problem, and that some form of early 
nonlocality may provide a resolution. The nonlocality could be generated by 
quantum nonequilibrium, or perhaps by some other means (other proposals 
include topological fluctuations [48] and an increased speed of light at early 
times [iS]). 

If we then assume — for whatever reason — that the pre-inflationary uni- 
verse was in a state of quantum nonequilibrium, the question is how the nonequi- 
librium will evolve in time, and in particular, whether any of it will survive until 
entry into the inflationary era. To address this question, let us first summarise 
what is known so far about relaxation in pilot-wave theory. 

As already mentioned, numerical simulations for simple two-dimensional sys- 
tems show an efficient relaxation, with an approximately exponential decay of 
the coarse-grained H- function H{t) [531 13S]- Specifically, for an ensemble of 
nonrelativistic particles in a two-dimensional box (on a static spacetime back- 
ground), with a wave function consisting of a superposition of the first 16 modes, 
it was found that H{t) « Hoe^'^^*" where, as discussed in rcf. [23j, the timescale 
tc coincides approximately with a theoretical relaxation timescale r defined by 
l/r^ = —{\/H)d?'H/dt^ [15], which under certain conditions may be roughly 
estimated as [T71[^l3n] 

r - At EE 1/Ai; , 

where A_E is the quantum energy spread and At is the usual quantum timescale 
over which the wave function evolves. 

Similar results have been obtained for nonrelativistic particles in a two- 
dimensional harmonic oscillator potential |39] , a case that has immediate impli- 
cations for the field theory of a single decoupled mode k. 

Writing 4* = V'kl'Zki, 91:2, i)>?, where x depends only on degrees of freedom 
for modes k' ^ k, equations (|lip and (|13|) imply that the wave function ip^^ 

^Note, however, that the existence of the puzzle depends on assuming a classical Friedmann 
expansion a oc f^/-^ all the way back to t = 0. 
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satisfies 



dt 2a3 \dql, dql 



while the de Broghe velocities for (gki, '7k2) are 



) ^k + ^afc' ('Zki + 9k2) ^k . 



(20) 



1 9sk . 1 9sk 
9ki — -^-^ , 9k2 — -^-^ (.^J-j 

(with i/ik = iV'kl s**'')- The marginal distribution ^^('Zki, '?k2, i) will evolve ac- 
cording to 

%+ E 7^(pAP^)-^- (22) 



As discussed elsewhere [25l [50] , these are identical to the equations of pilot- wave 
dynamics for an ensemble of nonrelativistic particles with time-dependent 'mass' 
m = , moving in the two-dimensional gki ^ 9k2 plane, in a harmonic oscillator 
potential of time-dependent angular frequency lj = k/ a. In the short-wavelength 
limit, Aphys << Artk • (where nk — n\^i + nk2 is the sum of the occupation 
numbers for modes kl and k2), and over timescales At = 1/Ai?k << 
(for which a is approximately constant), the above equations reduce to those 
for a decoupled mode k on Minkowski spacetime [10] . These limiting equations 
are in turn just those of pilot-wave dynamics for an ensemble of nonrelativistic 
particles of constant mass m = a'^ in a two-dimensional harmonic oscillator 
potential of constant angular frequency lo = k/a. The numerical results for this 
last case ^39, show that, in the Minkowski limit, for a decoupled mode k in a 
superposition of many different states of definite occupation number, one will 
obtain relaxation /3k('?ki, 'Zk2, i) |'0k(9ki, 9k2, 01^ (c"^ a coarse-grained level, 
assuming appropriate initial conditions), on a timescale Tk of order 

1 

Tk 



A^k 

If, in the Minkowski limit, relaxation occurs so efficiently for a single de- 
coupled mode, then we may reasonably expect that for a realistic entangled 
quantum state — in some pre-inflationary era — relaxation will occur at least 
as efficiently. One might then conclude that, even if there is initial nonequilib- 
rium, it will have relaxed away by the time inflation begins. However, before 
drawing definite conclusions, one must first consider the possible effect of spatial 
expansion on the relaxation process. One finds, in particular, that the character 
of the evolution can be very different in the long- wavelength limit. 

In the case of a decoupled mode on expanding space, described by the equa- 
tions ([^ - P^ . it is found [SU] that in the long-wavelength fimit, Aphys >> 
Arik ■ H~^, the wave function Tp^^ is approximately static — or 'frozen' — over 
timescales ~ H~^. Furthermore, one expects that the trajectories (qki {t), qk2{t)) 
will be frozen over timescales ^ H^^, in which case an arbitrary nonequilibrium 
distribution IV'kl^ will also be frozen over timescales ^ H^^. (This is 
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of course reminiscent of the freezing of super-Hubble modes in the theory of 
cosmological perturbations [Dll].) It then begins to appear possible that the 
normal process of relaxation to quantum equilibrium could be suppressed for 
long-wavelength modes in a pre-inflationary era, and that remnants of initial 
nonequilibrium could survive up to the beginning of inflation. 

That this is indeed possible has been shown jSU] by deriving a general and 
rigorous condition for the freezing of quantum nonequilibrium, a condition ap- 
plicable to an arbitrary time interval and to any (generally entangled) 
quantum state of a scalar field. (The condition may also be applied to mixed 
states and to interacting fields.) The condition is obtained by considering the 
displacements of the de Broglie-Bohm trajectories over the time interval [ti,tf]. 
It is found that, for a pure subensemble with (time-dependent) mean occupation 
numbers (fikr), nonequilibrium will be frozen (or at least partially frozen) for 
modes with wave number k if the time evolution of (nkr) satisfies the 'freezing 
inequality' [50] 




For a radiation-dominated expansion on [ti,tf], with a{t) = af{t/tfy/^, this 
inequality implies that (using (fikr) > 0) 

Aphys(i/) > 2TrHJ^ Hi flu) , (24) 

where Hj"^ = 2tf and the right-hand side is larger than Hj"^ if tf > {1.17)ti. 
Thus, in a radiation-dominated expansion, if the freezing inequality (j23p is sat- 
isfied, the corresponding modes must be super-Hubble [50] . 

We are now in a position to begin to address the question of whether or not 
very early nonequilibrium in a pre-infiationary era could survive until the onset 
of inflation itself. The above considerations show that, for short-wavelength 
modes, any initial quantum nonequilibrium is likely to be rapidly destroyed 
during a pre-inflationary phase. On the other hand, relaxation can be suppressed 
for long-wavelength modes — if the freezing inequality (|23p is satisfied — and 
it is then possible that these modes (for whatever fields may be present) will 
still be in nonequilibrium at the onset of inflation. 

Denoting for a moment the (approximately) constant Hubble radius during 
inflation by , relevant cosmological fluctuations originate from inside H^j . 
For some of these modes to be out of equilibrium, they must have evolved 
from modes that were outside the Hubble radius in the (presumably radiation- 
dominated) pre-infiationary phase. Therefore, for this scenario to work, some 
pre-inflationary nonequilibrium modes must enter the Hubble radius during the 
transition to the inflationary phase, and they must avoid complete relaxation to 
equilibrium by the time inflation begins. (As we shall see in section 6, relaxation 
does not occur during inflation itself.) Now, modes of physical wavelength 
Aphys = oA can enter the Hubble radius — a/d only if Aphys increases more 
slowly than does H~^, that is, only if the comoving Hubble radius H~^/a = 
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1/d increases — as occurs for a decelerating universe, a < (which, from the 
Friedmann equation a/a ~ — (47rG/3)(p + 3p), requires that the energy density 
p and pressure p satisfy p + 3p > 0) . 

During a decelerating pre-inflationary phase, then, any 'frozen' nonequilib- 
rium modes at super-Hubble radii can enter the Hubble radius. Once they do 
so, they are likely to begin to relax to equilibrium. For all modes that are inside 
at the onset of inflation, some time will necessarily have been spent in what 
might be crudely termed the 'relaxation zone', with Aphys ^ H~^, during the 
pre-inflationary phase. For example, for a radiation-dominated pre-inflationary 
phase (starting at some initial time ti) that makes an abrupt transition to an 
inflationary phase ait = tj, we have a = af{t/tfY/'^ and — 2t (on [ti,tf]), 
and a mode of comoving wavelength A enters the Hubble radius (aA ~ H^^) at 
a time tenter (A) ^ '^'f^'^/tfi so that the time spent in the relaxation zone is 

Airclax(A) = {tf — tenter) ^ tf{l — a'j / t'j) . 

There can be significant residual nonequilibrium at the beginning of inflation, 
provided the 'no relaxation' condition 

Atrelax(A) < t(A) (25) 

is satisfied, where t(A) is again a relaxation timescale as defined above. Because 
t(A) will depend on the wave functional, a proper calculation of r(A) requires 
a specific model of the pre-inflationary phase. 

Given a specific form for the function t(A), the condition (pS)) will deter- 
mine a range of wavelengths A for which residual nonequilibrium may reason- 
ably be expected to have survived from the pre-inflationary era. Because pre- 
inflationary modes with larger values of A enter the Hubble radius later and 
so spend less time in the 'relaxation zone', the condition ([25|) will presumably 
imply that residual nonequilibrium will be possible for A larger than some infra- 
red cutoff Ac. (The scenario might be improved if, during the transition from 
a pre-inflationary to an inflationary phase, the Hubble radius was a rapidly- 
increasing function of the scale factor {dH^^ /da >> H~^/a). For then super- 
Hubble nonequilibrium modes could be pushed far inside the Hubble radius in 
a short time.) 

We hope that future work, based on a speciflc pre-inflationary model, will 
yield a prediction for the infra-red cutoff Ac. It is of course possible that Ac 
will turn out to be so much larger than today's Hubble radius that it yields a 
negligible effect on CMB predictions. This remains to be seen. In this paper, 
the focus is on 'phenomenology': we simply assume that some modes could be in 
quantum nonequilibrium at the beginning of the inflationary phase, and we show 
how CMB data may be used to set experimental bounds on such nonequilibrium. 
Still, the above preliminary reasoning already suggests that if there is residual 
nonequilibrium from a pre-inflationary phase, then we should expect to find it 
at large wavelengths, beyond some cutoff Ac. 
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4.2 Possible Production of Nonequilibrium at the Planck 
Scale 

We have discussed whether nonequiUbrium might have survived into the in- 
flationary phase, on the assumption that there was nonequihbrium in some 
pre-inflationary era. Another question is whether nonequihbrium might be gen- 
erated during (or indeed even before) the inflationary era. 

The creation of quantum nonequihbrium from a prior equihbrium state is 
impossible in standard de Broglie-Bohm theory (leaving aside extremely rare 
fluctuations 1151), though it might occur in ahernative hidden- variables models 
— for example in models that deviate from quantum theory for processes taking 
place over very short timescales |51) . But even in de Broglie-Bohm theory, we 
do not know how to incorporate gravitation [37' . It is therefore conceivable that 
effects involving gravity are able to upset the equilibrium state. In particular, 
as has been discussed at length elsewhere, it is not unreasonable to propose that 
quantum nonequilibrium can be generated by the formation and evaporation of 
a black hole 22, 25j. 

This proposal is motivated by the (controversial) question of information 
loss in black holes. In the standard picture of black- hole formation and evapo- 
ration, it appears that a closed system can evolve from an initial pure state to 
a flnal mixed state, thereby violating ordinary quantum theory [52j . Further, 
because the final state describes thermal radiation that depends on the initial 
mass of the hole but not on the details of the initial state, it is impossible even 
in principle to retrodict the initial state from the final state. While Hawking's 
original argument for information loss remains controversial, a new approach to 
avoiding information loss invokes the possible existence of quantum nonequilib- 
rium in the outgoing radiation, which could then carry more information than 
can ordinary radiation in a conventional (mixed) quantum state. A mechanism 
for the creation of such nonequilibrium has been outlined [22j . involving an 
assumed nonequilibrium behind the horizon (presumably near the singularity) 
that is transferred to the exterior region by the entanglement between the in- 
going and outgoing modes of the Hawking radiation. A simple rule has been 
suggested, whereby the decreased 'hidden-variable entropy' (equal to minus 
the subquantum _ff-function ([2])) of the outgoing nonequilibrium radiation bal- 
ances the increase in von Neumann entropy ^vonN = — Tr(plnp) generated by 
the pure-to-mixed transition: 

A (^hv + ^voun) = . (26) 

Possible experimental tests of this proposal are discussed in refs. |22l [25] . 

It is sometimes suggested that, at the Planck scale, processes will occur 
involving the formation and evaporation of microscopic black holes. If one takes 
this (rather heuristic) picture seriously and combines it with the above proposal, 
one is led to the conclusion that quantum nonequilibrium will be generated at 
the Planck scale. During the inflationary phase, such processes might have 
an effective description in terms of nonequilibrium modes of the inflaton field 
at Planckian or trans-Planckian (physical) frequencies. One might reason as 
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follows. If a mode of comoving wavelength A once had a physical wavelength 
Aphys = aX < Ip near the beginning of inflation, one could assume that upon 
exiting the Planckian regime (that is, once Aphys becomes bigger than Ip) the 
mode will be out of equilibrium, having encountered some gravitational process 
that generates nonequilibrium while Aphys Ip- Whereas, modes that were 
never smaller than Ip will not encounter any such process. Roughly, one could 
model this by introducing a cutoff A[. such that nonequilibrium exists only for 
comoving wavelengths A < A^ {below the critical value A^, in contrast with the 
scenario in the preceding section). In addition to providing an estimate for A^, 
one also needs to estimate the degree of nonequilibrium, which for a given mode 
kr may be quantified by the relative (or hidden-variable) entropy 



An estimate for S'hv(fc) might arise from an application of (|26p in some form, 
though this remains to be studied. 

It is to be hoped that further development of this idea will lead to a detailed 
prediction for the form and magnitude of nonequilibrium for modes emerging 
from the Planckian regime. Pending such development, again, in this paper 
we restrict ourselves to using current data to set limits on any hypothetical 
quantum nonequilibrium that may be present during the inflationary phase. 

5 Measuring Primordial Quantum Fluctuations 

The above considerations suggest that, during inflation, some field modes may 
exhibit nonequilibrium fluctuations that violate quantum theory. Our aim in 
this paper is to show how to use CMB data to set bounds on such violations. 

We shall first recall how measurements of the CMB today allow us to infer 
statistical properties of inflaton fluctuations during the inflationary era. This 
involves working backwards from the CMB data, first to classical curvature 
perturbations in the early universe, and from these backwards even further to 
inflaton fluctuations during the inflationary phase. After having highlighted 
the key assumptions that are made in the standard treatment, we will be in 
a position to understand exactly how corrections to quantum theory during 
inflation are able to have an effect on the CMB. 

5.1 CMB Observations and Primordial Curvature Pertur- 



Employing angular coordinates {9, (j)) on the sky, CMB measurements provide 
us with a temperature function T{9, 0). Writing AT {9, (j)) = T{9, (j)) — T, where 
T is the average temperature over the sky, the temperature anisotropy may be 




(27) 



bations 
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decomposed into spherical harmonics 



oo 



T 



^ ^ atmYijn(e,(t)) 



(28) 



1=2 m=-l 



(where as usual we omit the dipole term) . A mode / corresponds to an angular 
scale « 6O7L 

A complete measurement of the microwave sky provides us with one function 
T(9, (p), or equivalently with one set {aim} of coefficients. In order to carry out 
a statistical analysis of {aim}, it is usually assumed (if only implicitly) that the 
observed T{6, 0) is a single realisation of a stochastic process, whose probability 
distribution P[T{9, 0)] (which may be thought of as representing a theoretical 
'ensemble of skies') satisfies the condition of statistical isotropy: 



for arbitrary angular displacements SO, Scj). This condition implies that, for a 
given I, each aim has the same (marginal) probability distribution pi{ai.m), with 
variance 



(the angular power spectrum, where (...) denotes an average over the theoretical 
ensemble) . 

Thus, given the assumption it follows that for each I we have what are 
in effect 21+1 independent realisations of the same random variable (with the 
same probability distribution). The observed quantity 



(constructed from measurements made on a single sky) then provides an unbi- 
ased estimate of the angular power spectrum Ci (that is, (c'^^^'^ = Ci), with a 
'cosmic variance' given by 



For large values of I, the quantity Cf ^ is an accurate estimate of Ci (that is, we 
expect to find Ci^^ ^ Ci). For small values of however, C^^^ is an inaccurate 
estimate of Ci . 

The observed CMB anisotropy is caused by classical inhomogeneities on the 
last scattering surface (when the CMB photons decoupled), and these inho- 
mogeneities in turn originate from classical perturbations that were present at 
much earlier times. In the long 'primordial' period between toxit(fc) and icntor(fc) 
(during which k « Ha, or Aphys >> H^^), the classical curvature perturbation 



p[T{e-5e,<p-5ct>)] = p[T{e,ct,)] 



(29) 




(30) 





(31) 




(32) 
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is time independent. (Here, ^^^i?k is the Fourier component of the spatial cur- 
vature scalar on comoving hypersurfaces, that is, on hypersurfaces with zero 
momentum density.) To a good first approximation we may ignore gravita- 
tional waves, in which case TZ^ is the only independent degree of freedom for 
the classical primordial perturbations. In terms of TZk, the aim may be expressed 
as [S3] 

a;™ = ^ / d^k T{k, 07^kl^™(k) , (33) 

where the transfer function T(fc, I) encodes the astrophysical processes that 
generate the temperature anisotropy. 

A given primordial curvature perturbation 7?.k (for all k) generates one set 
{dim} of temperature-anisotropy coefficients. A probability distribution P[7?.k] 
for IZk will generate a probability distribution P[{a/m}] for {aim}- If we make 
the assumption of statistical homogeneity, that P[7?.k] is translationally invariant 
— that is, in position space, P[7?.(x — d)] = P[}Z{x)] for arbitrary displacements 
d — it follows that (7^k7^k-) e~'^^+^'>-'^ = (7^k7^k-) and so 

(7^k7^^-) = 5kk•(|7^k|') . (34) 

From p3p and (|34p . the angular power spectrum (I30p may be written as 
where 

V^{k) ^ ^ (|7^k|') (36) 

is the primordial power spectrum. We shall assume, as is usually done, that 

l^kl^^ is a function of k only. 

Current measurements of the CMB show that V-R{k) w const, (an approxi- 
mately flat or scale- free spectrum) j54j . 



5.2 Inflationary Slow-RoU Predictions 

Standard inflation predicts an approximately flat primordial power spectrum 
Vn{k)- Let us briefly review how this comes about. 

An approximately homogeneous inflaton field (j)i^{t) -I- cj){x,t) (where is a 

small perturbation), with a potential V, has an energy density p w i^Q-f V^(0q). 
In the slow-roll approximation, p « V{(I)q) is approximately constant in time. 
The Friedmann equation (a/a)^ = (SvrG/S)/? then implies an approximate de 
Sitter expansion, a oc e^*, where H = ^ (87rG'/3)F((^o)- The time evolution of 
<^o is given by 

a dipr, 
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(where in the slow-roll approximation we may neglect the term ^q). The flatness 
conditions for V are e << 1, |?7| << 1 where 



1 f 1 dV 



) 



2 



1 1 (PV 

SttG V d(j)l 



e = 



WttG \ Vd(l), 



V = 



(37) 



The primordial perturbations are generated by quantum fluctuations dur- 
ing the slow roll. As a first approximation, the quantum fluctuations may be 
calculated for an eternal de Sitter expansion, and in this approximation one 
obtains an exactly scale-free primordial power spectrum. Corrections to this 
approximation yield small corrections to the scale- free result. 

The quantum theory of primordial perturbations has been developed in great 
detail TIIMIISS]. In the slow-roll limit {H 0), with V satisfying the flatness 
conditions, the inflaton perturbation = 4>{yi, t) evolves like a free massless field 
(until at least a few e-folds after toxit(fc) for the mode k). The quantised field 
(j) is usually assumed to be in the vacuum state. One may then use standard 
quantum field theory to calculate the probability distribution for the infiaton 
perturbation (j>. 

It is usually assumed that, a few Hubble times or e-folds after icxit(^) (that is, 
in the 'late-time limit'), the resulting quantum probability distribution for (j) may 
be regarded as a classical probability distribution over classical perturbations </>. 
This assumption has been justified by WKB-type classicality at late times [56] . 
by squeezing of the inflationary vacuum state [57|, I58j , and by environmental 
decoherence [521 HO]- The latter, in particular, seems to distinguish the field 
configuration basis (that is, the basis of eigenstates of the field operator (j)) as 
a robust pointer basis, where the relevant interactions are local in field space 
[Ml W5\ . The resulting distribution of field configurations is then, for practical 
purposes, indistinguishable from a classical distribution. Recent studies seem to 
confirm these conclusions: the pointer states consist (more precisely) of narrow 
Gaussians that approximate eigenstates of [5T] , and the locality of interactions 
in field space ensures that at late times the density matrix becomes essentially 
diagonal in the field configuration basis [62] • (See also ref. [63] for further 
discussion of WKB classicality in the late-time limit.) 

Given a classical inflaton perturbation 0, the corresponding curvature per- 
turbation is given by [1] 



where (fc) is a time a few e-folds after toxit(^)- The perturbation TZi^ is time- 
independent between i^(fc) and the approach to icntor(fc) (long after inflation 
ends), and is believed to seed what eventually grow into the dominant pertur- 
bations in the CMB. 

Note that the inflaton perturbation (j) is defined on a spatially flat slicing. 
(The inhomogeneous field (j) necessarily vanishes on comoving slices, since the 
momentum density —(j)V4> is by deflnition zero on such slices.) Then, in the 
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slow-roll limit H ^ 0, the effect of metric perturbations can be ignored [T]. The 
curvature perturbation TZ is defined on the comoving slicing. Thus, psp relates 
quantities defined on different slicings. 

The predicted (quantum-theoretical) primordial power spectrum, for 7?.k at 
t = t^, is then given by 



where 



(39) 



is the power spectrum of the infiaton fluctuations. 

As we have said, to a first approximation the infiaton fiuctuations are usually 
taken to be quantum vacuum fiuctuations in de Sitter spacetime. From the 
standard field operator expansion 
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(41) 

in terms of mode functions 

0+(x, « (^^ + e'('' "+'=/^'') (42) 

(solutions of ^ reducing to positive-frequency Minkowski modes in the short- 
wavelength fimit k/a >> H), the Bunch-Davies vacuum is defined by Ck |0) = 
(for all k). In this quantum state, the two-point (equal-time) correlation 
function is 

(oi0(x,t)^(x,t)io) = y: ^^^2^^"^' ^'"'^"^^ (43) 

k 

(where the first term in the numerator gives a Minkowskian contribution l/47r^a^|x- 
x1^). The quantum variance of each mode is given by the Fourier transform of 
the quantum two-point function, 

(Mqt-J^ J rf'xe-''-(O|,^(x + y)0(y)|O) , 

yielding 

<l*-l%T-5^|^(i + 7^) ^ («) 

The width decreases with time, tending to a finite constant. The power spectrum 
is 
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In the long- wavelength limit k/a « H (Aphys » H 
well outside the Hubble radius, we have 



471-2 



where the mode is 



(46) 



(If instead we set k = Ha, then V^^{k) = H^/2tt'^.) 

To a lowest-order approximation, then, the quantum fluctuations of the in- 
flaton field generate a scale-free spectrum of primordial curvature perturbations: 



1 



El 



(47) 
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These perturbations 7?.k c>c <j)y^ remain 'frozen' outside the Hubble radius until 
the time icntcr(fc) is approached. 

Because H and (jy^ are in fact slowly changing during the inflationary phase, 
higher-order corrections lead to a small dependence of V^{k) on k. 



6 Time Evolution of Nonequilibrium Vacua 

We now turn to the effect of quantum nonequilibrium on the predictions of 
inflationary cosmology. (A brief, preliminary account was given in ref. [25].) 
Our strategy is to consider nonequilibrium corrections to the lowest-order (scale- 
free) quantum spectrum, and then to compare these effects with the higher-order 
quantum corrections. 

In the Bunch-Davies vacuum, a mode kr has wave function ^^.j, = ^^.^((/kr, t) = 
\ipkr \ e*'''"' with a Gaussian amplitude 

^27rA 

of width 



IV'k.r = ^—e-''-/^^^ (48) 

k 



(contracting in time, and independent of r and of the direction of k) and with 
a phase 

---2imf^fc^ + ^W' (^«) 

where 

h(t) = - — - tan"^ — 

is independent of q^r- (It is readily verified that the above wave function 
V'kr('7kr, t) Satisfies the Schrodinger equation (|I6p for a mode kr, and that in the 
limit H — > 0, a ^ I one recovers the wave function ■(/'kr(9kr, i) c>c e^'^^^rg^^s'^* 
for the Minkowski vacuum.) 
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In the quantum vacuum, the q^r are independent random variables, each 
with a Gaussian distribution of zero mean. The width of each Gaussian decreases 
with time, approaching the asymptotic value H/ V2k^ (in the long- wavelength 
limit k/a << H). In the nonequilibrium (de Broglie-Bohm) vacuum, in con- 
trast, each Qkr evolves deterministically in time, and the probability distribution 
for each gtr depends on what the probability distribution was at some 'initial' 
time. 

The phase ([50|) implies a de Broglie velocity field 

dt ^ a3 dqur ^ H^a^ ' 

To solve ((5T|) for the trajectories qkrit), it is convenient to introduce the 
conformal time r/, defined by drj ~ dt/a. (For a cx e^* we have r/ = —1/Ha; as 
t runs from — oo to -|-oo, rj runs from — oo to 0.) In terms of rj, the equation of 
motion for gtr reads 

dll ^ 1 + k^TJ^ ' ^ ' 

which has the solution 



gkr(r/) = gkr(0)v/l + fcV • (53) 
The width of the packet is given by 

A^ = i^(l + fcV) . (54) 

An arbitrary distribution PkrC^kn ^) (generally 7^ iV'krl^kn '7)1^) necessarily 
satisfies the continuity equation 

5Pkr , 9 ( dq^r\ „ 



d-q dqi^r \ ^ drj 
which for the velocity field ([5^ has the solution 



Pkr('7kr, ??) = , ^ ^ Pkr(gkr/Vl + fc^'?^ 0) (55) 



VT+k-Tj 

for any given Pkr(9kr, 0). 

The time evolution amounts to a simple (homogeneous) contraction of both 
\ipkr\ ^'iid Pkr- Af times 77 < 0, IV'krl a contracting Gaussian packet of 
width Afc(77) ~ Ak{0)\/l + k'^rj'^, and in the late-time hmit 77 0, iV'krl^ 



approaches a static Gaussian of width Afc(O) = H/v2k^. At times 77 < 0, Pkr is 
a contracting arbitrary distribution of width D\^r{ii) = ^kr(0)\/l -I- k'^rf^ (with 
arbitrary Dkr(O)), and in the late-time limit 77 ^ 0, Pkr approaches a static 
packet of width Dkr(O) (the asymptotic packet differing from the earlier packet 
by a homogeneous rescaling of qkr, as in (|55p ). 
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For simplicity, we assume that (like A^) the nonequilibrium width Dkr is 
independent of r and of the direction of k, so that Dkr = Dk (t) . We then have 
the result 

— (const, in time) = Vf^fc) ■ (56) 

Note that, for each mode, the 'nonequilibrium factor' ^(fc) may be defined at 
any convenient fiducial time (in particular, not necessarily at the same time for 
every k). At least in this lowest-order approximation for the quantum state ^E*, 
it makes no difference whether we set the initial conditions for nonequilibrium 
at the same time for all values of k, or at different times for different values of 
k (for example at t{k) such that Aphys(fc) exceeds some critical value). 

7 Nonequilibrium Power Spectrum 

The above result for the nonequilibrium Bunch-Davies vacuum may be written 
as 

{\4>^.?)^W)^^m ■ (57) 

This implies that the nonequilibrium power spectrum for the inflaton fluctua- 
tions takes the form 

V^{k)^V'^^{k)ak) , (58) 

which for k/a « H reads 

^^(^) = ^^(^) • (59) 
The primordial power spectrum for the curvature perturbations is then 

Vn{k)=V^^{k)m , 



where r^^{k) is given by (|47p . Thus we have 



^0 J t.(fc) 



(60) 



In general, ^(fc) 7^ 1 and scale invariance is broken. In future work, along 
the lines outlined in section 4, we hope to be able to predict features of the 
function ^(fc). For the purposes of this paper, ^(fc) is (in principle) an arbitrary 
function to be constrained by observation. 

8 General Remarks 

Before considering how CMB data may be used to constrain the nonequilibrium 
function ^(fc), we make some general remarks on the above scenario. 
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8.1 Transfer of Microscopic Nonequilibrium to Cosmolog- 
ical Scales 

We saw in section 6 that, for each mode k during the inflationary phase, the 
respective widths Dk (t) and (t) of the noncquihbrium and equihbrium distri- 
butions remain in a fixed ratio Dk {t) / {t) — y^£,{k) over time. This holds in 
the approximation where the inflationary phase is treated as an exact de Sitter 
expansion. At least to a first approximation, then, we may conclude that quan- 
tum nonequilibrium (if it exists) will not relax during the inflationary phase, 
but is instead preserved over time. 

Furthermore, because of the exponential expansion of physical wavelengths 
Aphys during inflation, nonequilibrium (if there is any to start with) will not 
only be preserved but will also be transferred from microscopic to macroscopic 
scales. This 'magnification' of the nonequilibrium lengthscale is particularly 
striking in the late-time or large- wavelength limit Aphys >> H^^, where the 
de Broglie velocity field tends to zero for each mode, gkr ~* 0. In this limit, 
which takes effect a few e-foldings after the mode exits the Hubble radius, both 
p^^ and IV'krl^ become frozen. Once this happens, any difference between p^^ 
and iV'krI preserved, and is transferred to larger and larger lengthscales as 
the physical wavelength Aphys = a(i)(27r/fc) oc e^* of the mode gets larger and 
larger. The 'frozen' nonequilibrium then exists at a physical lengthscale that 
grows exponentially with time, from microscopic to macroscopic scales. 

Once inflation has ended, there will be a 'frozen' nonequilibrium distribu- 
tion of curvature perturbations TZ^ at macroscopic lengthscales. These per- 
turbations are then transferred to cosmological lengthscales by the subsequent 
(post-inflationary) Friedmann expansion. 

8.2 Quantum Measurement of the Inflaton Field 

As we saw in section 5.2, in the standard quantum theory of inflationary cosmol- 
ogy it is usual to assume that, during inflation, when the physical wavelength 
of a mode significantly exceeds the Hubble radius, the corresponding infiaton 
perturbation effectively 'becomes classical' — in the sense that the final quan- 
tum probability distribution for infiaton (and hence curvature) perturbations 
behaves to a good approximation like a classical probability distribution. As 
mentioned in section 5.2, various studies seem to confirm the validity of this 
assumption [56l ISll [Ml ESI [60l [Ml [62l [63] . In particular, the basis of eigenstates 
of the field operator ^ (suitably smeared with narrow Gaussians) seems to act 
as a robust pointer basis, so that the quantum distribution of field configura- 
tions is, for practical purposes, indistinguishable from a classical distribution 

[Miisniisi]- 

In the pilot-wave formulation of inflationary cosmology, there is a well- 
defined inflaton configuration or 'beable' (in Bell's terminology [64]) at all times, 
even before Hubble exit. In writing the formulas ((59]) and ((60]) . we have tac- 
itly identified the inflaton beable after Hubble exit with the 'classical' inflaton 
field after Hubble exit — where the latter generates the primordial curvature 
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perturbation via equation ([55]) . This identification merits some comment. 

Generally speaking, in pilot-wave theory, it is the precise value of the to- 
tal beable configuration that (together with the wave function) determines the 
outcome of a subsequent quantum measurement. However, as a rule, one must 
be cautious about identifying beable values with quantum measurement values: 
their relationship must be established on a case-by-case basis, through analysis 
of the particular measurement process that is taking place. As is well known 
PSI [M] , there are circumstances where a quantum measurement outcome does 
not provide a faithful record of the actual prior value of the beable (in which case 
the so-called quantum 'measurement' is in fact not a true measurement). For 
instance, in the pilot-wave theory of nonrelativistic particles, while the outcome 
of a quantum position measurement usually has the same value as the actual 
particle position prior to the measurement, for a quantum momentum mea- 
surement the outcome usually does not simply coincide with the prior particle 
momentum given by de Broglie's velocity formula. Instead, the quantum mo- 
mentum outcome depends on the initial particle position in a way that depends 
on the details of the measurement process. Thus, while quantum position mea- 
surements are usually 'faithful', quantum momentum measurements are usually 
not. 

Similarly, we expect that in the pilot-wave theory of fields, a quantum mea- 
surement of the field configuration will (usually) provide a faithful record of the 
value of the actual field beable appearing in the de Broglie-Bohm dynamics. For 
other measurements, however, this simple identification will not hold: instead, 
the outcomes will depend on the initial field beable in a way that depends on 
the details of the 'measurement' process. 

Now, in the case at hand, conventional analysis of the 'quantum-to-classical' 
transition during infiation indicates that the environment effects a quantum 
measurement of the inflaton field in the basis of field configurations [59l [60l [6T1 
162) . If this is correct, then we are indeed justified in our above identification of 
the de Broglie-Bohm inflaton field after Hubble exit with the 'classical' inflaton 
fleld after Hubble exit. 

Should the conventional analysis (for some reason) turn out to be incor- 
rect — in particular, if the 'quantum-to-classical' transition involves effective 
quantum measurements of the inflaton field in a basis different from the field 
configuration basis — then there will be a more complicated relationship be- 
tween the de Broglie-Bohm infiaton field and the emergent 'classical' infiaton 
field, a relationship that will depend on the details of the effective measurement 
process. There would then also be a more complicated relationship between 
the nonequilibrium distribution for the inflaton beable and the nonequilibrium 
distribution for the primordial curvature perturbations. 

8.3 Weak Dependence on Pilot- Wave Dynamics 

It is worth noting that the above results for the time evolution of nonequilib- 
rium vacua are only weakly dependent on the details of the de Broglie-Bohm 
dynamics. The results are in fact determined by just two features: (a) there is 
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a field beable (/)(x, t) whose time evolution is continuous and differentiable, and 
(b) the dynamics is 'separable', in the sense that for a product quantum state 
^[<^i ^] — n ■0kr('?kr, ^) the velocity of each component gkr is independent of the 

kr 

other QkrS. 

To see this, note that from (b) the evolution reduces to that of a collection of 
independent one-dimensional systems. Then, in each one-dimensional configu- 
ration space with coordinate gkr, the local conservation of quantum equilibrium 

dt dqi^r ' 

for some velocity field fkr = I'krl'Zkri t), uniquely fixes Wkr as 



1 d\i^Mr,t)\ 



VWr{qkr,t) ^ — — / rfg^^- 



2 



l'0kr(9kr,t)| Jq^^ dt 



(assuming that IV'krl^ vanishes at infinity), as follows immediately by integrating 
(|6ip with respect to the coordinate gkr, from some fixed value gkr to oo. 

Thus, for the case at hand, the assumption of a differentiable and separa- 
ble evolution fixes the de Broglie-Bohm velocity field uniquely. Note that this 
uniqueness arises only because the system reduces to a collection of independent 
one-dimensional systems. It is only in one dimension that the local conservation 
of quantum equilibrium fixes the velocity field. In two or more dimensions, other 
velocity fields are possible, distinct from that of de Broglie and Bohm |65j . 

The conditions (a) and (b) could certainly be violated in other hidden- 
variables theories. There might, for example, be no field beable (j){x,t) at all. 
Also, it is possible to have a pilot-wave-type theory with a non-separable dy- 
namics 66 . Still, property (a) might well emerge in some limit from a deeper 
hidden-variables theory. And property (b) seems desirable, even if not strictly 
necessary. In any case, our point here is to emphasise that (a) and (b) are the 
only features of pilot-wave dynamics that really enter into our considerations. 



9 Bound on Primordial Quantum Nonequilib- 
rium 

Let us now illustrate how the available data may be used to constrain the 
nonequilibrium function ^(k) appearing in the result (j60p for the primordial 
power spectrum, where the observed spectrum VTz{k) = 'P-^ {k)^{k) consists of 
the usual quantum contributions together with possible nonequilibrium correc- 
tions i^^l). 

It is currently a very active field of research to determine the fc-dependence of 
the observed spectrum V-ji^k), and to compare the results with the fc-dependence 
of the quantum-theoretical prediction (k). It is straightforward to reinter- 
pret these studies as effectively providing constraints on the nonequilibrium 
function ^(fc). 
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The observed spectrum Vn{k) is usually parameteriscd in terms of the spec- 
tral index n(fc), defined by 

and the running of the spectral index, n{k) = dn/dlnk. For n{k) approximately 
constant, it is convenient to write the power spectrum in the form 



Vn{k)^Vn{k^)[Yj ' ^^^^ 

where fco is some chosen reference or pivot point. (Note that the definitions 
(|62|) . (|63|) of n{k) generally agree at fc = fco only, and they agree for all fc if 
dn(k)/dk — 0.) The index n(fc) may be written as a Taylor expansion 

1 / fc 
n{k) = no + - In — 

2 \ko 

where no = n{ko) is the spectral index at fc = fco, and Uq = (dn/dlnk)^ is the 
running of the spectral index at k = kg. 

The observed values of n(fc), n'(fc) may be used to set bounds on early quan- 
tum nonequilibrium. To illustrate this, we shall consider a current best-fit value 
of no, 

no = OMOtrolt (64) 

at fco — 0.002 Mpc~^ |54| . (Adding nonequilibrium parameters would of course 
affect the best-fitting procedure, but the value ((64| suffices here for illustration. 
A best fitting of nonequilibriimr inflationary models to CMB data is outside the 
scope of this paper.) 

We have Vn{k) = "P^^ {k)^{k), where (k) is predicted by standard 
inflationary theory. One may adopt the following parameterisation: 



^r(fc)=^r(fco)(^^j , (65) 

where n^'^(fc) is the usual (quantum-theoretical) spectral index, and 



fc 



^(fc)-l 



m = C(fco) (^^ j , (66) 

where ^'(fc) is the 'nonequilibrium spectral index'. The observed index (minus 
1) is then a sum 

(n - 1) = (nQT - 1) + (,. - 1) (67) 

of contributions from quantum theory and from nonequilibrium corrections. 

In the exact limit iJ — > 0, we have nqx — 1 = 0; and in exact quantum 
equilibrium, we have — 1 = 0. Slow-roll inflation predicts a small tilt [Tl I53j 

nQ'^(fc) - 1 = -6e + 2r/ (68) 
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where, in the definitions (|37p of e and 77, the quantities V and dV/d(j)Q are 
evaluated at texit(fc) (for which k = aH). 

Defining i/q = v{ko) and n'^"^ = ?i'5"^(fco), we obtain a bound for |j/o — 1| on 
the assumption that \n^'^ — 1| is indeed significantly less than 1 (as predicted 
by inflation). Otherwise, in principle, both n^'^ — 1 and vq — 1 could be large 
— with comparable magnitudes and opposite signs — and the observed small 
value of their sum uq — 1 = — 0.04lg'Q|^3 could be an accident. We assume here 
that the observed small value |no — 1| ;5 0-1 is not due to such a 'conspiratorial' 
cancellation. Then, roughly, we may write (ag ain at ko = 0.002 Mpc~^) 

l«?^-l|<0.1, |z^o-l|<0.1. (69) 

The bound |i^o ~ 1| ^ 0.1 on the nonequilibrium index may be converted into 
a bound on the hidden- variable entropy S'hv(fc) — defined by ^7} — for modes 
with k close to fco = 0.002 Mpc"^. (As we have seen, S'hv(^) is the relative 
entropy of pj,.^ with respect to |'0krl i ^-^d is a natural measure of the difference 
between p^^ and IV'krl^-) We have ^{k) = I?^/A^, where Dk and are the 
widths of pj^j. and l^/^krl respectively. We know that IV'krI is a Gaussian packet, 
and that in the late-time limit = H^/2k^. For the purposes of illustration, 
let us model pj^.^ as a Gaussian (of width Dk). We then have 

5hv(fc) = i(l-e(fc)+lne(fc)) , (70) 

with ^(fc) parameterised by (I66p . If J^(fc) varies slowly, then close to ko we may 
write I'ik) ~ vq. Taking ^(fco) = 1 and assuming that [vq — 1| is small, we have 

^hv(fc)«-i(i^o-l)'ln'(fcAo) . 

Restricting ourselves to a range of k close to k^, such that |ln(fc/fco)| < 0(1), 
we then have 

\SY..m<\{v^-\f <\(r^ ■ (71) 

Note that approximate equilibrium in this region (close to fco) does not 
preclude large departures from equilibrium at much smaller or at much larger 
values of k. 

10 Possible Low Power Anomaly at Small / 

In the low-Z region (say I < 20), the angular power spectrum is dominated by 
the Sachs- Wolfe effect (resulting from non-uniformities in the local gravitational 
potential on the last scattering surface). 

In this region, T^(fc, I) = T:HQjf[2k/ Hq) [T], where is the Hubble constant 
today, so that (using p5|) ) 
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For Vnik) = const, we then have 



/' 

^0 



■oo 




1 



Ci cx 



2/(Z + l) ' 



so that l{l + 1)C; = const, at low / (the Sachs-Wolfe plateau), as seems to be 
approximately observed. 

It has been suggested that the data contain anomalously low power at small 
I, though this is controversial. If there is such low power, it could of coiuse 
be due to some inadequate processing of the data (such as in the modelling of 
foregrounds) or to some local astrophysical effect. Otherwise, the signal could 
be primordial in origin, reflecting an anomaly in the underlying spectrum V-iz{k) 
of curvature perturbations. In the latter case, the explanation might lie in some 
modification of the standard inflationary scenario, or in new physics. 

If there is a low-power signal at small I requiring new physics, then quantum 
nonequihbrium provides a possible candidate. Taking Vn{k) = {k)^{k), 
and assuming (to a first approximation) that 'P^^{k) = const., we may write 



If ^(fc) = 1 everywhere, then Ci/Ci = 1. A low power anomaly, Ci < C; , 
could be explained by having ^(fc) < 1 in some suitable region of fc-space. 
Because the integral is dominated by the scale k w IHq/2, a significant drop 
in Ci requires ^(fc) < 1 for fc in this region, that is, ^(fc) < 1 for wavelengths 
A « {4tt/1)Hq^ (comparable to today's Hubble radius). 

To have ^{k) < 1 for a primordial perturbation mode k means that the width 
Dk of the nonequihbrium distribution for the corresponding inflaton mode is less 
than the quantum equilibrium width A^. It is reasonable to expect this, if one 
accepts the scenario of section 2, according to which quantum noise arises from 
statistical relaxation processes (presumably taking place in the very early uni- 
verse). On this view, it is natural to assume that early nonequihbrium would 
have a less-than-quantum dispersion, or ^(fc) < 1 — as opposed to an early 
larger-than-quantum dispersion (^(fc) > 1) which, while possible in principle, 
seems less natural. Furthermore, we saw in section 4.1 that, in a supposed 
pre-inflationary era, relaxation to quantum equilibrium can be suppressed on 
large spatial scales, and one expects that at the onset of inflation nonequihb- 
rium is most likely to have survived at wavelengths A > Ac, where the value of 
Ac remains to be estimated (pending the development of an appropriate pre- 
inflationary model). Therefore, it appears that a dip ^(fc) < 1 in the power 
spectrum below some critical wave number kc = 2-k/\c could be naturally ex- 
plained in terms of early quantum nonequilibriiim surviving from a very early 
pre-inflationary era, though this possibility remains to be developed in detail. 

As for the possible production of nonequihbrium in the Planckian regime 
(section 4.2), in the absence of a more detailed model we arc imable to give any 
strong argument for ^(/c) < 1, as opposed to ^(fc) > 1, for modes with A < A^ 




(72) 



33 



(where remains to be estimated, see section 4.2). Again modelling 
Gaussian, the hidden- variable entropy 5hv(fc) for a single mode is given in terms 
of £_{k) by (|7(I|. For a given value of Shv{k) — perhaps set by some application 
of ([^ — equation (|7II)) possesses two solutions for ^(fc), one with ^ < 1 and one 
with ^ > 1. That is, the same nonequilibrium entropy can be achieved by both 
a less-than-quantum and a larger-than-quantum dispersion. On the other hand, 
from the behaviour of the function 1 — ^ + In^, one sees that the solution with 
^ < 1 always has a smaller value of I'C — 1| than does the solution with f > 1, 
that is, the solution with ^ < 1 has a dispersion that is closer to the quantum 
value. On this (weak) basis, one might suggest that f < 1 will be preferred. A 
stronger argument for ^ < 1 might come from a detailed understanding of the 
preservation of information by means of nonequilibrium noise suppression in the 
outgoing quantum state of an evaporating black hole. 

In any case, focusing here on the comparison with observation, let us consider 
the effect at low I of some simple examples of functions £,{k). 

As a first example, motivated by a possible long-wavelength suppression of 
relaxation at very early (pre-inflationary) times, we take S,{k) = for k < kc 
and ^(fc) = 1 for k > kc, where the simple cutoff is used to model a suppression 
of quantum noise at wavelengths A > Ac = 2t: /k^. We then have 

Ci-C?^ , [""dk 



QT 



f " rlk 

m + l) ^jf{2k/Ho) . (73) 



Again, the dominant scale is fc « lHo/2, and the correction to Ci will be sig- 
nificant only if the range of integration (0, fee) overlaps substantially with this 
scale — that is, fcc cannot be much smaller than IHq/2. 

Note that if, instead, we did take kc « IHq/2, the correction to C; would 
not only be small, it would be unobservable even in principle, because it would be 
smaller than the cosmic variance pip . For k « IHq/2 we have approximately 
jf{2k/Ho) « {2H\/i2l + ly.y {2k/Ho)^\ so that 

Ci-C^'^ . , ,^ f 2Hl \^ f2kcY^ 



This correction falls off rapidly with increasing Z, and is very small even for 
the lowest values of I: for example, even taking 2kc/HQ w 1, we find (C4 — 
C2'^)/C2'^ « —6 X 10^^. Because such corrections are much smaller than the 
cosmic variance ACf^^ /Ci = y^2/{2l + 1), they cannot be measured meaning- 
fully, even in principle. To obtain a measurable effect, the cutoff kc in (|73p must 
not be small compared to IHq/2. 

A second example is motivated by the possibility of gravitationally-induced 
nonequilibrium at small scales, at wavelengths A < A^.. If we assume that the 
nonequilibrium takes the form of noise suppression (^ < 1), one might model 
this again with a simple cutoff, taking ^(fc) = 1 for k < k'c and ^(fc) = for 
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k > k[., where k[. = 27r/A^. We then have 



Ci-C?^ , ,^ rdk 



QT 



m + j^^ -^jf{2k/H,) . (74) 



For a significant effect, the range of integration (fc^,oo) must again overlap 
substantially with the dominant region k « IHq/2 — which now implies that k'^ 
cannot be much larger than IHq/2. 

As a third example, we consider a power law 



.0-1 



m = ako) { ^ ) (75) 



(with constant index vq). From (j72p we then have 

2l{l + mko)(^) I dx jf{x)x^«-\ (76) 

where 



QT ^,^V^uj ^2fco 



dx jf{x)x''°- 



/¥r[(3-i^o)/2]r[? + (i.o-l)/2] 



/o 4 r[(4 - i.o)/2]r[? + (5 - i.o)/2] ■ 

Should the existence of a low-power anomaly be confirmed, one might try 
to match the anomaly with one of the above nonequilibrium spectra (I73p . (j74p 

or dze]). 

According to the analysis in ref. |67j . cutting off the power below a wave 
number fee ~ 3x 10~^ Mpc~^ (comparable to the inverse Hubble scale Hq — 2.4x 
10^** Mpc^^) slightly improves the fit to the three- year WMAP data, but the 
improvement does not seem large enough to justify any conclusion that such a 
cutoff really exists. Still, the possibility of reduced power at large scales is worth 
exploring, since it could originate from an early nonequilibrium suppression of 
quantum noise (as discussed in section 4.1). 



11 Non- Random Phases and Inter-Mode Corre- 
lations 

So far, we have considered only the angular power spectrum Ci of the microwave 
sky, and how this could be affected by nonequilibrium corrections to the pri- 
mordial (scalar) power spectrum V-jiik). Here, we consider how primordial 
non-Gaussianity could arise from early quantum nonequilibrium. 

The primordial curvature perturbations 7?,k are usually assumed to consti- 
tute a Gaussian random field, for which the power spectrum provides a complete 
characterisation of the statistical properties. The phases of Gaussian perturba- 
tions are randomly distributed, and there are no inter-mode correlations. 

In standard inflationary scenarios, the Gaussianity of 7?.k arises directly from 
the Gaussianity of the quantum vacuum fluctuations of the inflaton perturbation 
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(f)]^. (The Gaussianity of 7?.k is not, as is sometimes claimed, a mere consequence 
of the central limit theorem.) In the quantum Bunch-Davies vacuum, the infla- 
ton probability distribution at conformal time rj is given by 

kr 

where, as we saw in section 6, each ji/'krl^ ^ Gaussian of zero mean and 
width = {H^/2k^) {l + Pr]'^). The two-point function (O|0(x, 77)(/)(x, r;)|0) 
is given by The three-point function {0\(j){x,r])(}){x,r])(j){x,r])\0) vanishes, 

as do all odd-point functions. Higher n-point functions (for n even) reduce to 
sums of products of the two-point function, as expected for a Gaussian random 
field. In quantum equilibrium, then, the generation of primordial curvature 
perturbations TZu oc (j)^^ by inflaton perturbations is a Gaussian random process. 

However, as a general matter of principle, the primordial perturbations could 
be non-Gaussian. And if one considers quantum noncquilibrium for the inflaton 
field, there is no particular reason why the noncquilibrium inflaton fluctuations 
should be Gaussian. 

Wc have already seen that, in quantum noncquilibrium, the probability dis- 
tribution Pkrl'Zkr, ?7) for a single mode of the inflaton field need not take the 
quantum Gaussian form (j48p . Simple forms of non-Gaussianity include a non- 
zero skewness or kurtosis of Pkr(9kr,f?) (where the marginal Pkril'kr,v) for gkr 
may in general be obtained from a correlated joint distribution P[qkr, as dis- 
cussed further below). But non-Gaussianity can take on a wide variety of forms, 
and various measures of it have been proposed. Some workers have reported sig- 
nificant primordial non-Gaussianity in the CMB data |68| , while others maintain 
that the data are consistent with primordial Gaussianity [54) . 

Let us show how quantum noncquilibrium can result in non-random phases 
and inter-mode correlations for the primordial perturbations. 

The coefficients aim ~ \aim \ e'*^'™ in the spherical harmonic expansion (I28p 
are of course generally complex numbers, and their phases (pi^ contain a lot 
of information about the morphology of the temperature anisotropy AT(0, (j)) 
(see, for example, ref. [69]). Assuming again that the underlying 'ensemble 
of skies' is statistically rotationally invariant, the probability distribution for 
each ipij^ must be independent of m. For a fixed value of I, wc then have 
21 + 1 phases ipi^ with the same probability distribution pi{ipi„J, and for large 
I we may use the measured values of the ipim to probe (</?;„). At least to 
a first approximation, current data is consistent with pii^Pim) being uniform 
on the unit circle. According to the basic formula (I33p . each aim is a linear 
combination of all the curvature perturbation components 7?.k. And according 
to the inflationary result ([38]) . each TZit. is proportional to the late-time inflaton 
perturbation Thus, the phase iy9;„ of each aim is ultimately determined by 
the phases 6'k of all the inflaton perturbation components 0k = \4>-k \ c*^''- 

In quantum equilibrium, the inflaton phases 0k have a time-independent 
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distribution jO^"'"(0k) that is uniform on the unit circle: 

This follows immediately from ((48|) : the real and imaginary parts of (p-^ = 
(2^/2 (^ki + «9k2) have a joint Gaussian distribution a e~^''ki+9k2)/2'^fc that is 
always constant on circles centred on the origin in the complex ^^.-plane. 

In quantum nonequilibrium, the inflaton phases can at some initial (con- 
formal) time 7/j have an arbitrary distribution pij(0k,'7i)- Will the subsequent 
time evolution generate a late-time distribution that tends towards uniformity 
on the unit circle? Not in the approximation considered here. For the trajecto- 
ries Qkriv) = '?kr(0)-\/l + kP'rf' obtained in section 6 imply that 

6lk(??) - tan-i (gk2(r/)/gki(r/)) = tan^^ (gk2(0)/qki(0)) . 

Thus, during inflation, the phase of each inflaton mode is static, so that any 
initial nonequilibrium distribution (with non-random phases) will remain un- 
changed over time, p-^^{9]s_, rj) — p^^(9]s_, ryj for all values of conformal time 77. (In 
the complex ^^.-plane, the evolution of the joint probability distribution for gki, 
Qk2 amounts to a purely radial contraction with time, so that the distribution 
Pk(^k,?/) of phases is time- independent.) 

We conclude that the time evolution during the inflationary era does not 
scramble the phases of the inflaton perturbations. Any initial non-uniformity 
(or non-randomness) in the phase distribution will remain frozen, all the way 
to the late-time limit 77 ^ 0. It would be interesting, in future work, to explore 
how this could affect the phases iy9;„ of the measured coefficients aim in the 
temperature anisotropy. 

We now consider nonequilibrium inter-mode correlations. In section 6 we 
assumed, for simplicity, that the nonequilibrium distribution satisfied the fac- 
torisability condition (jlSp , so that the modes were uncorrelated even in nonequi- 
librium. However, in principle, correlations among modes are possible: in quan- 
tum nonequilibrium, the inflaton modes can be correlated even though |'I'|^ (for 
the Bunch-Davies vacuum) is a product. 

In terms of conformal time ?/, an arbitrary correlated joint distribution 
P[qkri v] "^ill evolve according to the continuity equation 

f^i:^(-^)-. 

Because the wave functional is still that of the Bunch-Davies vacuum, the ve- 
locity field dqkr/drj is still given by ()52p and the trajectories in configuration 
space are still given (mode by mode) by the result (l53|) . Given the trajectories, 
the general solution of ([77)1 may be constructed using the property that P/ l^fj^ 
is constant along trajectories (where this follows from the fact that P and l^fj^ 
obey the same continuity equation). Replacing the labels kr by a single index 
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n, we may equate 

P(<?i(0),g2(0),....,(?„(0),....,0) 



|^i(gi(0),0)|2 |V^2fe(0),0)^... K(g„(0),0)|2 
with 

Piqiiv),q2iv), ■■■■,qn{v), 



Using the trajectories qn {il) = 'Zn (0)^/1 + k^V^ and 

l^n(gn,r?)|' ^ 1 ^ A„(0) 

(where the width An{r]) is given by ([5^ ). we deduce that 

, p/^Ai(0) A2(0) A„(0) \^A„(0) 

This is an exact solution for the evolution of an arbitrary distribution, expressed 
in terms of the distribution P (qi, (72, Qn, 0) at conformal time 77 = 0. 

The possibility of nonequilibrium allows the distribution P[q]ir,T]i] some 
initial time ?7j to be, in principle, anything at all. To narrow down the range 
of possibilities, one might impose the requirement of statistical homogeneity, 
P[^(x — d),r7j] = P[0(x),77j] (for arbitrary spatial displacements d). 

Clearly, allowing non-random phases and inter-mode correlations in the in- 
flationary vacuum opens up a number of novel possibilities. Indeed, the subject 
of non-Gaussianity for quantum nonequilibrium states deserves to be developed 
in more detail. 



12 Conclusion 

We have shown how inflationary cosmology (assuming it to be essentially cor- 
rect) may be used to test the validity of quantum theory at very short distances 
and at very early times. In particular, we have considered the possible effects 
of quantum nonequilibrium, as described by the hidden-variables theory of de 
Broglie and Bohm, during the inflationary phase. We have shown, by means of 
simple examples, how CMB data may be used to set bounds on nonequilibrium 
deviations from quantum theory. 

As for the possible origin of such deviations, we have outlined a scenario 
where quantum nonequilibrium during the inflationary phase arises from relax- 
ation suppression (for long-wavelength modes) in a pre-inflationary era. This 
scenario suggests that primordial nonequilibrium could set in above some infra- 
red cutoff Ac (though the value of Ac remains to be estimated). We have also 
considered the more speculative possibility that nonequilibrium could be gen- 
erated during the inflationary era, by novel gravitational effects at the Planck 
scale. 
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We have for the most part discussed quantum nonequiHbrium corrections 
to the primordial (scalar) power spectrum V-jiik)- A preHminary discussion 
was also given showing how primordial non-Gaussianity (in particular non- 
random phases and inter-mode correlations) could also arise from early quantum 
nonequilibrium . 

One important topic we have not discussed is how quantum nonequilibrium 
would affect the consistency relation between the power spectra for scalar and 
tensor perturbations. This relation is especially interesting because it relates 
fluctuations for distinct degrees of freedom, and because it is independent of 
the form of the inflaton potential. Presumably, quantum nonequilibrium would 
in general have different effects on different degrees of freedom, resulting in a 
violation of the consistency relation. 

If standard inflationary cosmology is essentially correct, then observations 
of the CMB have already confirmed — to a first approximation — the validity 
of the quantum (Born-rule) prediction for the inflaton power spectrum during 
the inflationary phase. More accurate measurements of the primordial power 
spectrum will enable us to set unprecedented bounds on violations of quantum 
theory, at very short distances and at very early times. And close scrutiny of 
other possible features, such as various forms of non-Gaussianity, will provide 
further tests of basic quantum predictions. 

Should inflation be very firmly established, and should it be found that the 
predictions of quantum theory continue to hold good at all accessible length- 
scales during the inflationary era, then this would constitute considerable evi- 
dence against the hypothesis of quantum nonequilibrium at the big bang (though 
of course, nonequilibrium from an earlier era might have simply not survived 
into the inflationary phase). Furthermore, it would rather undermine the view 
that quantum theory is merely an effective description of an equilibrium state. 
In principle, one could still believe that hidden variables exist, and that the 
hidden-variables distribution is restricted to quantum equilibrium even at the 
shortest distances and earliest times. But in the complete absence of nonequilib- 
rium, the detailed behaviour of the hidden variables (such as the precise form of 
the trajectories in de Broglie-Bohm theory) would be forever untestable. While 
exact equilibrium always and everywhere may constitute a logically possible 
world, from a general scientific point of view it seems unacceptable, and the 
complete ruling out of quantum nonequilibrium by experiment would suggest 
that hidden-variables theories should be abandoned. 

On the other hand, a positive detection of quantum nonequilibrium phenom- 
ena in the early universe (or indeed elsewhere 25J) would be of fundamental 
interest, opening up a new and deeper level of nature to experimental investi- 
gation. 
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